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BACKGROUND OF THE INVENTION 
Field of the Invention 
The present invention relates to a method of generating a hybrid grid 
allowing modelling of a heterogeneous formation crossed by one or more pipes. 
5 The method of more particularly applied to formation of a grid suited to an 

underground reservoir crossed by one or more wells, in order to model 
displacements of fluids such as hydrocarbons. 




Description of the Prior Art 

1 0 Grid generation is a crucial element for the new generation of reservoir 

simulators. Grids allow description of the geometry of the geologic structure 
studied by means of a representation in discrete elements wherein simulation is 
performed according to a suitable numerical pattern. Better comprehension of 
physical phenomena requires 3D simulation of the multiphase flows in 

1 5 increasingly complex geologic structures, in the vicinity of types of singularities 
such as stratifications, faults, pinchouts, channels and complex wells. All this 
complexity has to be taken into account first by the grid which has to reproduce 
as accurately as possible to geologic information in its heterogeneous nature. 
Grid modelling has made great advances during the past few years in 

20 other fields such as aeronautics, combustion in engines, structure mechanics, 

etc. However, the gridding techniques used in the other fields cannot be applied, 
as they are, to the petroleum sphere because the professional constraints are not 
the same. For example, in reservoir simulation, the numerical patterns are 



2 

constructed from control volumes in order to better respect the mass 
conservation in the case of transport equations of hyperbolic nature. The grid 
must be a "block-centered" type grid, that is the nodes must be inside each 
layer and the boundaries of each block must follow the interface between the 
5 layers. Now, if this constraint was not taken into account, the nodes would 

naturally be placed along the faults and along the stratification boundaries. The 
consequence of this would be that these interfaces would pass through the 
control volume that is used. The saturation, constant in the control volume, could 
not consider discontinuity and the results would not be accurate. It is therefore 

10 necessary to develop new techniques that are better suited to requirements of 
the petroleum field. 

Cartesian grids, which are commonly used in current commercial 
simulators, are not suitable for solving these new problems posed by the 
development of petroleum reservoirs. Cartesian grids, based on parallelepipedic 

1 5 elements, do not allow representation of such complex geometries. 

There is a well-known method of generating structured 3D hexahedral 
grids of CPG (Corner-Point-Geometry) type which considers the geometry of the 
bodies, It is described in French patent 2,747,490 (US patent 5,884,564) filed by 
the assignee and also in the following publication: 

20 - Bennis Ch. Et al. "One More Step in Gocad Stratigraphic Grid 

Generation: Taking into Account Faults and Pinchouts" ; SPE 35526, Stavanger, 
1996. 
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This grid type is more flexible than Cartesian grids because it consists of 
any hexahedral elements that can be degenerated. It strictly respects the 
horizons, the faults and it allows representation of certain nonconformities such 
as pinchouts because its construction is based on these elements. However, this 
5 type of grid does not allow solution of all the geometric complexities such as, for 
example, circular radial grids around complex wells. It is possible to form 
separately the grid of the reservoir and the girds around the wells but it is difficult 
to represent several objects in the same CPG type reservoir grid because of 
connection problems linked with the structured nature of the grid. 

10 Another approach is also known where 3D grids only based on tetrahedral 

Delaunay elements, with a circular radial refinement around the wells, are 
automatically generated. The advantages of such an approach is that it is 
entirely automatic and does not practically not require the user's attention. 
However, this method has drawbacks which make the results obtained difficult to 

15 use: 

there are on average five times as many grid cells as in a CPG type grid 
for the same structure, which is very disadvantageous for simulation calculations, 

unlike the structured grids which are easy to visualize, to explore from the 
inside and to locally modify interactively, it is very difficult and sometimes 
20 impossible to properly control the tetrahedral grids because of their size and 
especially because of their non-structured nature. This poses problems for 
validating the grid from a geometric point of view as well as for understanding 
and validating the result of a simulation of this type of grid. 
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Other approaches are also well-known, which allow generation of grids, 
notably grids based on control volumes generated from a triangulation (VoronoT 
and CVFE), associated with techniques of aggregation of the triangles (or 
tetrahedrons) into quadrangles allowing the number of grid cells to be reduced. 
5 Although promising results were obtained with these new grids, precise 
representation of the geologic complexity of reservoirs and wells remains a 
subject for research and development. In fact, these grids are rather 2.5D (i.e. 
vertically projected) and their 3D extension appears to be very difficult. Despite 
their hybrid aspect, they remain entirely unstructured and would therefore be very 
10 difficult to manage and to handle in real 3D. Furthermore, taking account of real 
3D faults and deviated wells would greatly increase this difficulty. 



SUMMARY OF THE INVENTION 
The method according to the invention allows generation of a hybrid grid 
15 applicable to a heterogeneous formation crossed by at least one pipe of known 
geometry (such as an underground reservoir crossed by one or more wells), in 
order to form a model representative of fluids flows in the formation in 
accordance with a defined numerical pattern, the structure of the formation being 
known a priori from available data acquired by means of in-situ measurements, 
20 analyses and/or interpretations of formation images (seismic images for example, 
in the case of a reservoir). 

The method comprises associating a first structured grid for gridding the 
formation with consideration of discontinuities thereof with second structured, 
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radial type grids for gridding of the zones around the wells, these second grids 
allowing consideration of constraints linked with flows in the wells, and non- 
structured transition grids between the first grid associated with the formation and 
the second grids associated with the wells. 
5 Gridding of the heterogeneous medium is obtained, for example, by 

importing each second structured grid into a cavity formed in the first structured 
grid, the size of the cavity being sufficient to allow formation of a non-structured 
transition grid between the first structured grid associated with the formation and 
the second structured grid associated with each well. 

10 The non-structured transition grid based on any polyhedrons and 

canonical polyhedrons such as pentahedrons, tetrahedrons, pyramids, etc., can 
be formed by consideration of constraints linked with the numerical pattern. 

The non-structured transition grids are advantageously modelled with the 
structured well grids by applying a technique known in the art, referred to as the 

1 5 "generalized map" technique, the grid of the formation being structured 
metrically, globally or in faulted blocks. 

The global hybrid grid is thus obtained by combination of several grid 
types : a structured reservoir grid, a radial grid around each well, also structured, 
and non-structured transition grids which connect the previous two grid types. 

20 Each one of these grids has its own formation, representation and exploration 

methods. The structured aspect is thus degraded only at the points where this is 
strictly necessary. This "object" approach affords both the advantage of 
structured grids for control and comprehension of the reservoir and the flexibility 
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of non-structured grids in complex zones. Complexity is introduced only where it 
is strictly necessary. The independence of these gridding modes therefore 
allows separate extraction, management and representation of the well grids and 
the interstitial grids includes in the reservoir grids. 
5 Using a reservoir simulator of a well-known type, such as ATHOS™ or 

SCORE™ for example, for a reservoir provided with a hybrid grid obtained by 
means of the method, allows production simulations to be performed. 



BRIEF DESCRIPTION OF THE DRAWINGS 
10 Other features and advantages of the method according to the invention 

will be clear from reading the description hereafter of non limitative examples, 
with reference to the accompanying drawings wherein: 

Fig. 1 shows a diagrammatic example of a hybrid grid of a reservoir 
crossed by two wells, consisting of a first structured grid for the reservoir, a 
15 second structured grid for the zones around the wells and transition grids 
between the two grid types; 

Fig. 2 shows an example of a structured grid of a faulted reservoir; 
Fig. 3 shows an example of a radial grid around a vertical well; 
Fig. 4 shows an example of grid of a horizontal well; 
20 - Fig. 5 shows, in 2.5D, an example of a gridded reservoir where cavities 
are provided for gridded wells, before the stage of creation of non-structured 
interstitial grids intended to connect them together; 
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Fig. 6 shows five wells provided each with a radial grid, included in a 
gridded reservoir, by means of non-structured transition grids based on any 
polyhedral grid cells; 

Fig. 7 shows separately the structured reservoir grid with the cavities 
5 provided therein in order to include the additional elements : gridded wells and 
interstitial grids; 

Figs. 8A, 8B, 8C, 8D show elementary well grids illustrated individually 
according to different modes; 

Figs. 9A, 9B, 9C, 9D show various elementary transition grids, illustrated 
10 individually according to different modes, also, allowing their integration into the 
reservoir grid; 

Fig. 10 shows an example of a hybrid reservoir grid with transition grids 
consisting for example of pentahedrons, between the reservoir grid and several 
well grids; 

15 - Fig. 1 1 shows a model for a matrical representation of a structured grid 
around a well; 

Fig. 12 is a graphic representation of a connection between strands, within 
the scope of the modeling technique referred to as the "generalized map" 
technique used for generating non-structured grids; 
20 - Figs. 13A, 13B are graphic representations of connections by means of 
simple arcs; and 

Figs. 14A, 14B are graphic representations of connections by means of 
double or triple arcs respectively. 



DETAILED DESCRIPTION OF THE INVENTION 
Modelling of the reservoir is obtained by combining elementary grids of 
different types. Each grid is considered to be a full object with its own data 
5 model, its own generation methods and its own representation methods. 
Generation is carried out in stages with addition/subtraction of grids. 

1 ) In order to represent the reservoir as a whole, an i, j, k structured grid 
of a known type , referred to as CPG (Corner Point Geometry), as described in 
the aforementioned French patent 2,747,490 for example is used . The reservoir 

10 can be faulted with downcreep of a block in relation to the other. The major 

horizons and faults are first modelled by continuous surfaces from data resulting 
from an interpretation of seismic images of the subsoil or from data obtained 
during drilling (well markers). The geologic structure is then divided into faulted 
blocks resting on these surfaces. These blocks are individually gridded, then 

1 5 reassembled. Gridding of a block first consists in gridding the edge surfaces, 
then the inside is populated by transfinite interpolation of the edge surface grids. 
Relaxation techniques are then applied to the edge surfaces and to the inside so 
as to harmonize and the regulate the grid. The grid thus obtained strictly 
considers the horizons, the faults and allows representation of nonconformities 

20 such as pinchouts. It meets all the constraints of geologic nature. 

2) A well trajectory is drawn synthetically or imported. A structured radial 
grid is then generated around each well in order to take into account particular 
constraints linked with the flows in the vicinity of the wells. 
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In the example shown in Fig. 3, the structure grid around a vertical well is 
of circular radial type. It is also a CPG type grid. Its generation first consists in 
sampling a disc at r, 9 in the horizontal plane. The 2D grid thus obtained is then 
projected vertically upon the various layers of the reservoir grid. Here, the i, j, k 
5 of the matrical structure correspond to the samplings at r, G and z respectively 
(see Fig. 11). 

The grid around a horizontal well (Fig. 4) is i, j, k structured, it is of the 
same type as that of the reservoir, except that a well cannot be faulted. It is also 
obtained by projecting vertically upon the various layers of the reservoir grid a 2D 
10 grid belonging to a horizontal plane. 

3) This radial grid is the inserted around around each well in the global 
reservoir grid. A cavity is therefore first created in the reservoir grid by 
deactivating all the grid cells in contact with well grid cells (Figs. 5, 6). The space 
freed between the reservoir grid and the well grid must be sufficient to allow 

1 5 convenient formation of a transition grid. It can represent for example about the 
equivalent of two grid cell layers. 

4) A non-structured transition grid is then generated in this cavity (Fig. 7) 
in order to connect the structured radial grid around the well to that of the 
reservoir best respecting the constraints linked with the numerical pattern. The 

20 user can deactivate the grid of a well any time by reactivating the grid cells of the 
corresponding cavity in the reservoir grid. 

The transition grid can for example consist of polyhedrons with any 
number of sides or canonical polyhedrons (tetrahedrons, pentehedrons, 
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pyramids, etc.) according to the numerical pattern used, without the overall 
hybrid approach proposed being affected. 

Example of modelling of a hybrid grid 

5 The reservoir grid and each well grid are modelled, for each faulted block 

of the reservoir, by matrical structures of points or cells each comprising eight 
points. Because of the structured nature of the grids, the topological links 
between the various grid cells are implicitly contained in the matrical structure. 
Transition grids are more difficult to manage because of their non- 
10 structured nature and because they can contain polyhedral grid cells whose 
number of sides varies from one cell to the other. 

An advantageous solution for facilitating management of this new grid 
type, allowing to browse it and to surf it efficiently, uses the topological model 
referred to as the "generalized maps or G-maps". This model known and for 
1 5 example described by: 

Edmont J.: "A Combinatorial Representation for Polyhedral Surfaces", 
Notice Amer. Math. Soc, 7, 1960, or by: 

Fortune S., 1992: Voronoi Diagrams and DelaunayTriangulations, pp. 225- 
265 of D.Z. Du & F.K. Hwang (eds.), Computing in Euclidean Geometry, 2 nd edn. 
20 Lecture Notes Series on Computing, vol. 4, Singapore, World Scientific. 

Generalized maps are based on a formal algebraic approach that is briefly 
described hereafter. 
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In 3D, the elements which constitute a generalized map are (D, a 0 , ch, a 2 , 
a 3 ), where D is a finite set of elements called strands, and elements {a} are 
involution on D type functions, associating the strands two by two at most, which 
are therefore conveniently referred to as links. Figs. 1 1 to 14 show specific 
5 geometric representation examples. Link a 0 is in the form of a dotted segment 
(Fig. 11) and links Ch, a 2 and 03 in the form of arcs, respectively simple (Fig. 12), 
double (Fig. 13) and triple (Fig. 14). 

According to another known approach, generalized maps are considered 
as graphs whose strands from the nodes and the links form the arcs : link ct 0 
10 between two strands can be used for representing the edge of a side, links ai for 
connecting two edges of a side, links a 2 for linking two sides of a cell together 
and links 03 for connecting two cells together. 

The model of generalized maps involves a small number of formal objects 
and an operation which, by associating additional information with a topology, 
1 5 allows location of the objects defined in space and accounting for their 

appearance, which is referred to as "plunge" and, in the present case, "plunge in 
a 3D space". 

It affords the advantage of being independent of the dimension of the 
objects. All the objects can be represented with the same data structure and 
20 handled with the same methods. This approach makes it possible to handle 
objects created with heterogeneous topological models. It is therefore well- 
suited for implementing the method according to the invention with creation of a 
non-structured grid linking two different structured grids together. 
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The generalized map technique for modeling the transition grid is applied 
by creating a number of objects of different types which refer to one another. 
These objects define the topological network and its various plunges in a 3D 
space. Specifically, in order to allow browsing the grid, a topological network is 
5 constructed parallel to the geometric data commonly handled in a grid, the 
points, the edges, the sides and the cells. Furthermore, crossed links are 
established between the topological network and geometric data. 

Objects 

The various objects handles within the scope of the application performed 
here of the generalized map method are as follows: 

1) The Transition Grid object which contains all the topology, the 
geometry and the physical data. It is a GMap type object which represents the 
topological network and of a Plunge type object which defines the plunge in the 
physical world according to the application. 

2) The GMap object 

The topological model is entirely contained in a graph of a list of Strands 
connected to one another. Any operation performed on the generalized map 
20 amounts to an operation of the Strands network. The GMap object type has 
methods allowing easy circulation in the topological network representing the 
grid, that is to go from one Strand to another. 
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15 



3) The Strand object 

Each Strand is defined by four references to other Strands (corresponding 
to links a 0 , ai and 03) and by four other references to the plunge in the 3D space, 
notably at a Point, an Edge, a Side and at a Cell to which it is connected. 

5 

4) The Plunge object 

It is defined by four lists : 

a list of Points (the grid points), it is the plunge of dimension 0 of the 

GMap, 

10- a list of Edges (plunge of dimension 1 of the GMap), 

a list of Sides (plunge of dimension 2 of the GMap), and 
a list of Cells (plunge of dimension 3 of the GMap). 

The Plunge object also contains its own methods of creating and handling 
the data it contains according to the use that is made thereof. The GMap is 
15 created from its Plunge. 



5) The Point object 

A Point is defined by its coordinates x, y, z and by a list of attributes, 
notably scalar or real petrophysical values that are associated therewith. 

20 

6) The Edge object 

It is defined by a reference in the GMap to a Strand which represents an 
end of the Edge. This gives a preferential access to the topological network and 
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simultaneously allows going from the plunge to the strands graph. For example, 
link a 0 of the Strand in question leas to the Strand representing the other end of 
the Edge. It is furthermore defined by a list of attributes, notably scalar or real 
petrophysical values that are associated therewith. 

7) The Side object 

This type of object allows handling directly the interfaces between the cells 
as well as the outer sides of the grid. A side is defined by a reference in the 
GMap to a Strand which represents a vertex of the edge polygon of the Side. 
This also gives a preferential access to the topological network. The Strands 
representing the other vertices of the polygon are accessible by iterative 
applications of the relation a 0 oai by starting from the initial Strand and eventually 
coming back to this same Strand. It is also defined by a list of attributes relative 
to the Side (for example, scalar or vector physical values). 

8) The Cell object 

The type of Cell object is defined by a list of references in the GMap to 
Strands, each on representing a half-Side of the Cell. This allows access to the 
topological network from the Cells. It is also defined by a site (coordinates of the 
center of mass of the Cell) and by a list of attributes specific to the Cell or to its 
site (scalar or vector petrophysical values for example). 

Graphical representation and exploration 
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A graphical representation is a very efficient and even essential means for 
controlling and validating the construction of a grid and the simulation results. 
Concerning the constructions, the geometry of the grid generated is generally 
first visually controlled. If this is not sufficient, local or global quality criteria with 
5 which statistics are established can be calculated and visualized on the grid by 
means of a color scale. Flow simulation calculates the variations with time of 
certain petrophysical parameters by accounting for the hypotheses that initially 
condition the flows. Simulation validation also involves visualization of these 
parameters on the grid (preferably by means of a colour scale). As the involved 
10 grids are 3D grids, tools allowing exploration of the grid from the inside by visual 
browsing are required. Graphical representation and browsing in the grids, 
presented hereafter, are a good illustration of the flexibility and the modularity of 
the hybrid approach proposed and of the efficiency and the adequacy of the data 
model selected. 

15 The hybrid grid is considered as a set of independent entities. The 

elementary grids is constantly seen in a main window. The user can select at 
any time an elementary grid and see it with its specific methods in a secondary 
window which contains only the selected elementary grid. Actions on the 
elementary grid have automatically repercussions on the entire hybrid grid seen 

20 in the main window. An elementary grid can thus be seen and explored as a full 
entity and it can be viewed in the global context. The visualization methods differ 
according to whether the elementary grid is structured (reservoir grid and well 
grids) or not (transition grids). 
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Examples of functionalities specific to structured grids 

In the case of a structured grid (reservoir and well), visualization is simple 

and conventional. It consists in two main functionalities: 
5 - visualization of the external envelope of the grid with the possibility of 

peeling it in the 3 directions i, j, j separately, 

simultaneous or separate viewing of three matrical cell slices i=cste, j=cste 

and k=cste, with the possibility of moving slices in the block. 



10 Examples of functionalities specific to non-structured grids 

In the case of non-structured transition grids, other, more elaborate 
visualization modes are preferably used. Five functionalities are mainly used: 

viewing of the external envelope with the possibility of concentric peeling, 
topologically speaking, 
15 - viewing of the cells crossed by cutting plane orthogonal to an axis of 
coordinates x, y, z or any axis, 

viewing of the trace of the cells on the cutting plane, 
viewing of the grid sites when they are initially given and 
viewing of the grid cells in full or scattered mode. 
20 Of course, for the two grid types, it is possible to view a property or scalar 

value by means of a color scale. 

All these functionalities require easy and optimum viewing of a non- 
structured grid. This is possible by using the generalized map technique. 
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Figs. 7 to 10 clearly illustrate the potential afforded by the hybrid grid that 
is harmonious integration of a structured grid following a topological model (well 
grid) into another structured grid (reservoir grid) following a different topological 
model, by means of a non-structured transition grid. The independence of these 
5 models therefore allows extraction and separate representation of the well grids 
and of the interstitial grids included in the reservoir grid in order to represent, 
handle and explore this type of data. 



